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ABSTRACT 

Aims. We investigate the electromagnetic interaction of a relativistic stellar wind with a planet or a smaller body in orbit around the star. This 
may be relevant to objects orbiting a pulsar, such as PSR B1257+12 and PSR B1620-26 that are expected to hold a planetary system, or to 
pulsars with suspected asteroids or comets. 



Methods. We extend the theory of Alfven wings to relativistic winds. 



(Q Results. When the wind is relativistic albeit slower than the total Alfven speed, a system of electric currents carried by a stationary 
$-H , Alfvenic structure is driven by the planet or by its surroundings. For an Earth-like planet around a "standard" one second pulsar, the associated 
CZ3 ' current can reach the same magnitude as the Goldreich- Julian current that powers the pulsar's magnetosphere. 
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1. Introduction 

Precise pulsar timing measurents proved that the pulsars PSR 
B1257+12 and PSR B1620-26, host planets, at d i stance s 
of order of an astro nomical unit (IWolszczan & Frail I (Il992b . 
iThorsett et al. (1993)). Moreover, accretion discs are expected 
to form at some phase of the evolution of neutron stars in a bi- 
nary system, possibly giving birth to second generation planets. 
Small bodies, such as planetoids, asteroids or comets may also 
orbit pulsars and occasionally fall on them. 

Circum-pulsar objects move in the centrifugally driven rel- 
ativistic pulsar wind. The angular velocity Q* of a rotating neu- 
tron star typicall y is larger than 10 rad.s" 1 . The star behaves 
like an antenna ( Deutschll 19551) emitting by magnetic dipole 
radiation a power which causes it to lose rotational energy at a 
rate 

E rot = -M/Q»f2„ = 47T 2 M / / , /P 3 , 



(1) 

where Mj x (2/5)MR 2 is the neutron star's moment of inertia 
dLvne & Graham-Smith||l998l) . This power can be compared 
to the orbital energy of a pulsar planet. Consider, for exam- 
ple, the case of PSR B 1257+ 12 and its planet "a". Using the 
data in Tables Q] and [2 we can estimate the moment of iner- 
tia to be Mi = 1.x 10 38 kg.m 2 , the star's rotational energy 
loss to be E ror = -2. x 10 27 W, and the planet's orbital en- 
ergy to be E G = GM t M P /2a = 4. x 10 32 J, where a stands 



for the planet's semi-major axis . The planet and its environ- 
ment intercept a fraction larger than or of order of jiRpliAna 2 ) 
of this power. This captured power is used to heat the planet, 
to generate the current system described below and to work 
on t he planet's motion as desc ribed in the accompanying pa- 
per dMottez & Hevvaertsl201ll) . hereafter (MH2). The planet's 
radius can be infered from its mass by assuming a terrestrial 
density (5000 kg rrT 3 ). The intercepted power E capt is : 



-'capt 



(Rp \ 2 
~2a\ 



(2) 



For planet "a" of PSR B1257+12, we find, E cap , = 1.74 10 18 
Watts. If a substantial part of this power goes into performing 
work on the planet's motion, an orbital evolution time scale of 
about 8. 10 6 years may be expected. This is a short time scale 
by astronomical standards. It scales with the mass Mp of the 

1/3 

planet as M p ~ . 

In this paper, we examine the interaction of these planets 
with the magnetized wind of their pulsar and discuss in (MH2) 
the effect of this interaction on the long term evolution of their 
orbital elements. 

The energy flux carried by the wind of an ordinary star, 
such as the Sun, is small enough to have but negligible effect 
on the orbits of its planets. Stellar winds usually are asymptot- 
ically super-Alfvenic, by which we mean that they eventually 
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become faster than the total Alfven speed, associated with the 
modulus of the magnetic field. As a result, the planets are pro- 
tected from a direct contact with the wind by a bow shock. 

Pulsar winds are much different. In a first approximation, 
oblique rotator pulsars may be regarded as magnetic dipoles ro- 
tating at high angular velocity in vacuo, which causes a low fre- 
quency and large amplitude electromagnetic wave to be emit- 
ted, the wave character of which reveals itself in the wave 
zone, beyond the light cylinder, of radius r = c/Q*. In this 
zone, magnetic field lines become spiral-shaped and the az- 
imuthal component of the magnetic field decreases as ~ r' 1 
while its radial component decreases as B r ~ r~ 2 . The mag- 
netic field becomes mostly azimuthal at large distances and 
at the equatorial latitudes were the planets are expected to be 
found. Aligned rotators essentially are rotating unipolar induc- 
tors which generate, beyond the light cylinder, an highly rela- 
tivistic MHD wind in which the magnetic field also becomes 
predominantly azimuthal at large distances from the rotation 
axis, the field components and B, varying with distance es- 
sentially as indicated above, at least near the equatorial plane 
and when the flow is close to being radial. The power emit- 
ted by such objects is also of order of E rot , given by Eq. (Q]i, 
because the Poynting flux emitted through the light cylinder 
is comparable for wind or wave emission. We only consider 
aligned rotators in this paper. Their analysis is simplified by 
the fact that the magnetic field of the wind observed in the 
planet's frame is close to being time-independent. The origin, 
acceleration and structure of pulsar winds are not fully under- 
stood. A number of models have been proposed in the literature 
dMichellll969t iHenriksen & Ravburn|ll97U Icontopoulos et al 



1999; Michel 2005; B ucciantini et al.l2006h . In spite of their di- 
versity, they all converge on the fact that the wind is dominated 
by the Poynting flux, although observations indicate that they 
turn into matter-d ominated high energy flows at large distances 
dKirketal., 2009). At distances of order of an astronomical unit, 
pulsar winds are expected to be still Poynting-flux-dominated. 
This means that the electromagnetic energy density ~ Bq/jUo, 
fio being the magnetic permeability of vacuum, is much larger 
than the plasma energy density yopoc 2 , yo being its Lorentz 
factor and po the rest mass density of this supposedly cold 
wind. These quantities refer to the observer's frame. In such 
circumstances, Alfvenic perturbations propagate at a phase ve- 
locity close to the speed of light (equation (0). Therefore, in 
spite of being highly relativistic, the wind flow may neverthe- 
less be sub-Alfvenic, in the sense defined above. A planet in 
a Poynting-flux-dominated wind may remain unscreened from 
the wind by a bow shock and thus enter in direct contact with 
it. 

The interaction of a planet with a sub-Alfvenic plasma flow 
has been considered for moderately magnetized non-relativistic 
flows in connection with the interaction of the satellite Io with 
the plasma and magnetic field present in Jupiter's magneto- 
spheric environment. This interaction is driven by the induc- 
tive electromotive field which results from the motion of the 
satellite across Jupiter's corotational magnetic field and plasma 
flow. The satel l ite act s as a (uniformly moving) unipolar indue - 



explored case of such a structure in the universe up to now, 
and it has bee n the object of recent studies conce rning its over- 
all structure dChust et alJl2005l iHess et al.ll2010t), the p ossibil- 
ity of particle acceleration (IHess et al.ll2007bl 12 009b) a nd its 



consequences on the radio emissions dOueinnec & Zarkal 19981 



Hess et al. I l2007al 12009a). In the present paper, we develop a 
theory that generalizes some of Neubauer's results to the case 
of highly magnetized (B^ » /iopoyoc 2 ) and relativistic plasma 
flows with Lorentz factors yo >> 1> when a planet immersed 
in the magnetized pulsar wind acts as a unipolar inductor and 
generates two stationary Alfvenic structures that emerge from 
the planet and extend far in the wind. 

2. A unipolar inductor in the pulsar wind 

Let us consider a planet orbiting a pulsar in the relativistic flow 
of the emitted wind. Different reference frames can be involved 
in the description of the fluid motion. The frame where the neu- 
tron star is at rest is the observer's frame. We denote it by Ro- 
The planet velocity being small compared to the wind veloc- 
ity, we may consider the planet to be at rest with respect to 
the neutron star, except when discussing the planet's motion. 
The reference frame Ro can then also be regarded as being the 
planet's rest frame. Quantities observed in this frame are de- 
noted by letters without any superscript, such as p or v. The 
unperturbed wind's instantaneous rest frame in the vicinity of 
the planet is the "wind's frame" Rw- Quantities observed in this 
frame are denoted by letters with a prime superscript, such as 
p'. An index refers to quantities associated with the unper- 
turbed wind. The unperturbed wind velocity in Ro is Vo, its 
associated Lorentz factor is jq, the unperturbed magnetic field 
is Bo and the wind's density is po. The wind plasma being sup- 
posedly cold, this mass density reduces to the (apparent) rest 
mass density of particles. In the wind frame Rw, the unper- 
turbed wind density is its proper rest mass density p' Q and its 
magnetic field is B Q . They both differ from po and Bo. In the 
presence of a perturbation, the instantaneous rest frame of the 
fluid is not Rw but another rest frame, Rf . Quantities observed 
in this instantaneous rest frame are indicated by a subscript F. 

The perturbation generated by the planet in this flow is 
time-dependent in the wind's frame. Since in the tenuous and 
highly magnetized pulsar wind the formal Alfven velocity Ca 
may exceed the speed of light, the derivation of the propaga- 
tion velocity V' A of Alfvenic perturbations must take into ac- 
count the displacement current. In the wind's rest frame : 



V. 



'-2 



-2 



'-2 



c 2 + f*op' Q /B' 



o ■ 



(3) 



tor. 



Neubauer (119801) derived a nonlinear theory of this inter- 



action. The Alfven wing connecting Io and Jupiter is the only 



If the flow is faster than the Alfven velocity Va, the planet is 
preceded by a shock wave that defines a confined area where 
the flow is strongly modified. But if the flow is sub-Alfvenic, 
the planet is directly in contact with the wind. In the instanta- 
neous rest frame of the wind plasma we expect, following the 
MHD approximation, the electric field to vanish. In the planet's 
rest frame an electromotive field E = -vxB is generated, where 
B and v are the magnetic field and fluid velocity in the planet's 
frame. 

Using a model pulsar wind, it is possible to make an es- 
timate of the electric field. One of the simplest such models 
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was developed by Michel lll969h . It includes the magnetic field 
and the star's rotation, wind particles masses, but neglects the 
dipole inclination and the gravitation and assumes the wind 
flow to be radial. The unperturbed magnetic field Bo only has 
a radial poloidal component B' and an azimuthal toroidal com- 
ponent fig. The neglect of the latitudinal component B^ — 
is justified near the e quatorial plane and far e nough from the 
inner magnetosphere (Conto poulos et al. I ll999h . For bodies or- 
biting near this plane, as circumpulsar planets probably do, 
this geometrical restriction is unimportant. Other models for 
stead y state axisymmetric winds have been de veloped since 
then JBeskin et aflll998l: iBucciantini et alj|2006l) . The station- 
ary equations of an axisymmetric perfect MHD flow admit a 
set of integrals of motion along stream lines, such as the mass 
flux / and the magnetic flux *P, which, for cold radial flows and 
radial poloidal fields, are defined by: 



/ = yop' v r r 
*P = r 2 B' n , 



(4) 
(5) 



where all quantities refer to the unperturbed wind as seen in the 
observer's frame. The MHD app roximation E q + To X Bo = 
and the Faraday equation imply (Meste ll96ll) : 



Bl = K 



.v -Q,r 



(6) 



In order to judge whether the planet is in superalfvenic motion 
with respect to the wind, we should assess whether the mod- 
ulus of its velocity in the wind's frame, which equals that of 
the wind in the observer's frame, Vo, is faster or not than the 
propagation speed in the wind's frame, V' A , defined by Eq. (3). 
The square of the ratio of these two velocities is given, for an 
asymptotically radial wind, by : 

v2 ( 

(7) 



The wind being supposedly radial, vt vanishes. This simplifies 
equation (6) which also shows that at distances much larger 
than the light cylinder radius, the radial component of the mag- 
netic field can be neglected compared to the azimuthal one. The 
Alfven speed c A calculated in the observer's frame is given, us- 
ing equations (4) and (5) to express B' {) and p°, by: 



(8) 



The square of the Alfven speed c' A in the wind's frame is a 
factor yo smaller than the square of the Alfven speed c A in the 
observer's frame because po = yop' and B^ = yoB'^. The 
latter relation is a result of the vanishing of the electric field in 
the wind's frame and of the velocity vo being perpendicular to 
the azimuthal field. Then: 



K 2 = ["i 



1 + 2*1^ 

cr \ c 



where <tq is the magnetization parameter: 

Q 2< P 2 



cr = 



Mo/c 



3 - 



(9) 



(10) 



In highly relativistic Poynting-flux-dominated outflows, as is 
expected for pulsar winds, cr » 1. The models show that 

1/3 

asymptotically the Lorentz factor yo approaches cr ( j and thus 
(vo/c) 2 approaches (1 - o~ 2 ^ 3 ). From equation (0, it is found 
that, to lowest order in an expansion in cr ', the asymptotic 
value of M' A is (1 - l/(2o\^ 3 )). Therefore the pulsar wind re- 
mains slower than the total Alfven speed (3). 

Since this result applies at distances much larger than the 
light cylinder radius, it is essentially valid wherever planets 
may be found orbiting. We therefore consider that the plan- 
ets detected around pulsars orbit in a sub-Alfvenic relativistic 
pulsar wind. 

In the numerical computations, we are less subtle and state 
that at the planetary distances, vq ~ c and Va ~ c, every time it 
makes sense. Then, in first approximation, the unipolar induc- 
tor electric potential drop U (along the 9 axis, i.e. perpendicular 
to the orbital plane) is 



U = 2R P E = 2R P E e = 2R P v r B* 



2R P Q*\> 



(11) 



The Fig. [T] provides a representation of the geometric configu- 
ration of the unipolar inductor. 

The rotation period of the pulsar PSR 1257+12 is P — 
0,006s, which corresponds to O* = 1010,397 rad.s 1 and 
its surface magnet ic field is estimated to B t 8.8 x 10 8 G 
(Tay lor et al.ll2000i) . We assume a star radius R t — 10 km and a 
mass M* = 1.4M (see Table [TJ. Data concerning this planet 
can be read in table|2] We assume an Earth-like density of 5000 
kg.m 3 . The flux Y = R 2 t B* = 8.8 x 10 12 Wb. The semi-major 
axis, the planetary radius Rp = (iMplAnpp) 1 ^ and an estimate 
of U are given for each planet "a", "b", "c" orbiting this pul- 
sar in Tables [2] and [3] It can be seen that the inductor electric 
potential drop U (from pole to pole along the planet) is of the 
order of 10 12 V. 

For the pulsar PSR 1620-26, the rotation period is P — 
0,011s, Q* = 567 rad.s -1 , the surface magne tic field is esti- 
mated to B t = 3. x 10 9 G dTavlor et al.1 
that R* 



2000l) . Still assuming 
3 X 10 13 Wb. This pulsar 



10 km, we find that ¥ 
has a white dwarf companion star. The neutron star mass can 
1.35M Q dThorsett & Arzoumanianlll999 



be estimated to M, 



Sigurdsson et al ]l2003h . The planet is more distant from its star 
than in the case of PSR 1257+12 but the planetary radius is 
larger. The resulting potential drop U is still of the same order 
of magnitude. 

We show in sections[3}j5]that a planet in the wind generates 
two current systems that propagate far in space, forming a so- 
called Alfven wing. 

3. Equations of special-relativistic ideal MHD 

In this section, the equations up to Eq. ( 1231 are general to spe- 
cial relativity, and valid for quantities defined in any inertial 
reference frame. We use notations without prime and subscript. 
After Eq. (l23l , quantities without prime and subscript refer, as 
in the previous sections, only to the observer's frame Ro de- 
fined at the begining of section [2] 

A condition for MHD to be valid is that the typical length 
scale L relevant to the flow be much larger than the particles 
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Larmor radii p F = yv ± /(jj c , where v± is the typical velocity 
perpendicular to the magnetic field in the fluid's frame, a> c the 
gyrofrequency and y the particle's Lorentz factor in the same 
frame. The energy of particles in a pulsar wind is a priori very 
high since the wind's bulk Lorentz factor yo may be as large 
as 10 s to 10 7 . If however, in the wind's frame, a significant 
part of this energy resided in perpendicular motions, particles 
would very quickly loose it by synchrotron radiation, even in a 
moderate magnetic field. Therefore, we may consider the pul- 
sar wind particles to have negligible Larmor radii in the rest 
frame of the wind's bulk flow, so that MHD is applicable on 
almost any scale. 

The fourth component of position in space time is x 4 = 
ct, t being the time measured in seconds in the given refer- 
ence frame and c the speed of light. Greek indices label four- 
dimensional coordinates and components and latin indices la- 
bel Euclidean three-dimensional ones. The metric tensor of 
Minkowskian space is diagonal, with components 7744 = + 1 and 
711 = 1I22 = t?33 = -L Vfi designates the partial derivative with 
respect to the space-time coordinate xf. The notation V desig- 
nates the three-dimensional nabla operator. We use the dummy 
index rule. Special-relativistic MHD equations associate fluid 
equations with Maxwell's equations, in which the displacement 
current and Poisson's equation should be retained. The fluid 
equations consist of a conservation equation for particle num- 
ber, valid in the absence of reactions among particle species, 
and of the four components of the energy and momentum con- 
servation equations The law of conservation of particle number 
is written as: 

^n F uT) = 0. (12) 

The density n F is the proper spatial number density of particles, 
that is, their density measured in the instantaneous rest frame 
of the fluid. The four components u M are those of the dimen- 
sionless four-velocity of the fluid in the considered rest frame: 



Poynting vector divided by c and its space-space components 
form a second rank tensor of the three-dimensional Euclidean 
space, the Maxwell stress tensor: 



2 2/l/o 



2 ^ = = BB 

6 - e EE 



(16) 



The symbol 6 represents the second rank unit tensor. The mat- 
ter energy-momentum tensor of a cold pressureless fluid can be 
written, in the absence of non-ideal effects such as viscosity, as: 



(17) 



In the absence of internal heat, the number conservation equa- 
tion (fT2l reduces to a conservation equation for proper mass, 
since in this case pf = mrif, m being the rest mass of par- 
ticles. Equation (fT2l can then be written in three-dimensional 
notations as: 

dyp F 



dt 



+ div (yp F v) = . 



(18) 



Here, we need not solve the energy conservation equation be- 
cause the medium is regarded as cold. We only consider the 
spatial components of the equivalent equations (TPfl i or (HBl . 
From equation (TBI we get, denoting the charge and current 
density by p e and j respectively: 

— (y 2 p F v) + div (y 2 p F vv) = p e E + j x B . (19) 
From equation (fT31 l we get the equivalent equation: 



E xB 

— \y*p F v + - 

dt \ p Q c z 



div [y 2 p F vv + Mj = . 



(20) 



The components of the Maxwell stress tensor M may be trans- 
formed to account for the perfect MHD relation 



(yv/c, y) . 



(13) 



E + vxB =0. 



(21) 



The Lorentz factor y refers here to the bulk fluid motion. The 
number density in the observer's frame is the fourth component 
of the density-flux four- vector, n F y. The equations of conser- 
vation of energy and momentum of matter are lumped in the 
four-tensorial equation: 



v fi 1 m Jet 



(14) 



The T^'s are the components of the energy-momentum tensor 
of matter and the components f v em are those of the electromag- 
netic force density four-vector. By using Maxwell's equations, 
this four- vector can be written in conservative form and equa- 
tion (TBI) can be given the form: 



V7T + O = o. 



(15) 



The r^'s are the components of the electromagnetic energy- 
momentum tensor. They can be expressed in terms of the elec- 
tromagnetic field strength tensor F or in terms of the electric 
and magnetic fields observed in the chosen reference frame. 
Its time-time component is the electromagnetic energy den- 
sity, its space-time components are the three components of the 



The components of the electric field and those of the Maxwell 
stress tensor are easily calculated in a frame where the x-axis 
is taken to be along the direction of the fluid velocity v. Some 
algebra then yields the following expression, where the indices 
t (transverse) and / (longitudinal) refer to a component of a 
vector perpendicular or parallel to the velocity v of the fluid: 



B] - 1 Bf = BB B 2 vv v 2 B,B, 



M=— l -6+-r — -8 + — — + — 



2p y 2 2p p Q p Q c 



c 2 po 



(22) 



The magnetic field evolution equation, deduced from Faraday's 
equation and the perfect MHD relation, writes: 



dB 

— = curl(v x B) . 
ot 



(23) 



The non-relativistic theory of lNeubauer ( 1980 ) considers only 
the Alfvenic wake of the satellite. Indeed, fast MHD distur- 
bances propagate isotropically in the low-Bp limit and decrease 
in amplitude with distance from the source (Bp is the ratio of the 
plasma pressure to the magnetic pressure). Such disturbances 
do not create any concentrated current system. The slow mode 
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propagates in a low-ftp plasma much slower than Alfenic per- 
turbations and the associated disturbances, though channelled 
by the magnetic field, soon become spatially separated from the 
Alfenic wake. Actually, slow mode perturbations barely prop- 
agate at all in a cold pulsar wind and carry negligible current. 
The effects of compressive perturb ations in the non-relativistic 
situation have been discussed by Wrig ht & Schwartzl (119901) 
who have shown that compressive plasma wave modes, though 
necessarily excited, contribute one order of magnitude less to 
the current flow and to the energy budget than shear Alfven per- 
turbations do. We therefore follow Neubauer in concentrating 
on purely Alfvenic motions. We assume that the fluid motions 
triggered when the planet passes by are sub-relativistic in the 
rest frame Rw of the unperturbed wind and that the changes of 
the formal Alfven velocity are similarly sub-relativistic. We do 
not assume however that the formal Alfven speed calculated 
from the total field, B' /(fiop') 1 ^ 2 , is less than the speed of light. 
Our assumption of non-relativistic Alfvenic motions in Rw im- 
plies that div v' = 0. It does not follow however that the same 
relation also holds true in the observer's frame Ro because the 
wind flow is relativistic in this frame. In Rw, the electric terms 
of the Maxwell stress tensor are by of order v' 2 /c 2 less than 
the magnetic terms and can be neglected. Accounting for the 
relation div v' = 0, equation d23l can be written as: 



— -(B'-V')v' =0. 

df 



(24) 



Considering Eqs. (fJTJ and d22l . the equation of motion ( 1201 
becomes, in Rw, 

d / , , fl'V (v'-B')B' 
— 'P v +— T 



df 



p. a c^ p. a c<- 
a/B' (v'-B')v'' 



= 0. 



(25) 



4. An Alfvenic first integral 

The computation of the non relativistic Alfven wing given by 
Neubauer (1980) is based upon the fact that, in simple non- 
linear Alfvenic wave motions, the velocity V, = v - s Ca is 
a first integral. In this relation, v is the fluid's velocity, Ca the 
vectorial Alfven velocity associated with the perturbed field, 
and the sign s = +1 depends on the sense of propagation of the 
perturbation. This relation can be transposed in a differential 
form as c/B = kd\, where k = (yuop) 1 ^ 2 . Moreover, the modulus 
of the magnetic field is a time invariant and, because of the 
uniform boundary conditions B = Bq, this modulus is constant 
over the whole space. 

As Neubauer, we look for a solution where c/B' = kd\' and 
where B' is constant. Then, setting b' = B' /B' and 



2x1/2 



A = (fi p' + B' z /c z ) 
we can write the Eqs. (l24l and (l25l in the form 



df 



B' 2 (b' ■ V')b' 



B' 



dVM 



df 



- + (b'-V')(v'v') 



^ - (b' • W = 0, 

df 



(26) 

■■ (27) 
(28) 



where vj = v' • b'. The two first terms of Eq. (l27l are analo- 
gous to those of classical MHD, when Va « c. The two others 
are specific to fast variations of the electric field, when the dis- 
placement current is taken into account. We solve this system 
in the following way: we first ignore the last two terms of Eq. 
(fJTJi and solve the resulting system (I27l-(l2"8ll. which can then 
be written as: 



dv' 

A 2 fi' 2 (b'-V')b' =0, 

df 

dW 



df 



(W • V')v' = 0. 



(29) 
(30) 



We then check whether the last two terms in Eq. (127} vanish or 
are negligible. These terms are: 



B' 2 



dv'„b' 



df 



+ (W ■ V'Xv'vO 



(31) 



Considering Eqs. d29l and (1301 . the proportionality of c/B' and 
d\' is obtained for k = sA, and these two equations become 
equivalent to 



B'db' = sAd\' 

db' , , , , 
sA— - B'(b' ■ V')b' = 0. 

df 

Eq. (l32l has a first integral V, 

v: = v' - s 



B' 



(piop' +B' 2 /c 2 ) 1 ' 2 ' 



(32) 
(33) 

(34) 



In the region not perturbed by the planet, v' = 0, and B' = W Q , 
Therefore, 



' A' 
v^^fi'-b'-B^). 



(35) 
(36) 
(37) 



In order of magnitude, A ~ B' /c, and the approximation V « 
c implies that ||B' - By|| << B' Q . Therefore, this is equivalent 
to an hypothesis of linea r perturbat i on. Th is remark also holds 
for the computations of uSfeubauerl (119801) . At this stage, it is 
necessary to evaluate the terms in Eq. (|3~TV Using the above 
relations, these terms are 



c \ dt c 



o ' 



(38) 



where bg is the direction of the unperturbed magnetic field in 
the wind's frame of reference. The quantity in Eq. (l38l does not 
vanish but it is negligible, as long as v' << c, in comparison to 
the terms in Eq. (|29l , which are of the order of (B' 2 /c)db'/dt'. 
Therefore, in the linear approximation, the solution given by 
Eqs. (132H37t is correct. 

The first integral of Eq. (l34l must now be computed in the 
observer's frame of reference Ro- It will be used in the next 
section to derive the value of the electric current associated to 
the Alfven wings. Let us write the velocity and the magnetic 
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field as the sum of a longitudinal component, parallel to v<), and 
a transverse vector component. 



B, = 



B, 



VqV 


B, 


vg 




(- 


VoVo 


6- 




V 


-I 



B. 



(39) 
(40) 



The transforms of the velocity and of the magnetic field are 
, v/ - v + v,/y 



1 - v • Vo/c 2 

b; = b,, 

D , L /, v-v \ B-vo \ 



(41) 
(42) 
(43) 



where jq is the Lorentz factor of the unperturbed wind. For the 
transform of the magnetic field, Eq. (l2"TT i has been taken into 
account. The Eq. d34l i becomes 



v. 



v/ - v + v,/y 



1-m 



-iro 



B ■ V » I T> 



(44) 



For a first order development, we note Vi the velocity perturba- 
tion and Bi the magnetic perturbation, 



V = Vo + Vi B = Bo + Bi . 

To the first order, Eq. (l44l becomes 

V. = yl{vu+— )- ^(Ba + Bid 

\ 70 / A 

vo • vi B • v 



(45) 



A 



( 1 



7o 



■ (Bo/ + Bi,) - B 



■ Vi, 



(46) 



which can be rewritten 

V' s / s B • v \ Xn s s 

-T =Vi;+ T^B,- — 



Considering the purely geometric relation 

s B • v s VpVo 

Ac 1 Ac 1 



(\_ 



Vi • v 



B,. 



(47) 



Then we can define a new first integral vector, 



v - - B 

a 



c 2 c 2 



= v + 



Bq-Vq 



where the constant number a is defined by 

. B ■ vo 
a = A — s — . 



, (49) 



(50) 



In spite of a greater complexity, this vector presents some anal- 
ogy with the first integral found in the non relativistic case by 
Neubauer (written in the first lines of the present section). This 
will be used in the next section to derive the current flowing 
along the Alfven wing, in a similar way to those developed in 
Neubaueri (1 19801) . The parameter A, given in Eq. d26b . can be 
expressed as a function of the unperturbed wind parameters, 



A = [fi p' + c- 2 Bl + c- 2 y 2 Bl] 



1/2 



(51) 



Since V • v' =0, the density p' is invariant along any line of 
flow; that is why we have noted it p' a instead. We set 



B 2 /fiop' , 



Al 



(52) 
(53) 
(54) 



These are mixtures of the magnetic field in the frame Ro of the 
observer and of the density in the proper frame Rw of the wind 
in terms of which the factor A can be written as: 



-4 



HoPo 



1 + 



2 2 
C AlT 



-At 



y 2 c 2 



1/2 



(55) 



5. Current carried by an Alfven wing 

Neubauer (1980) has shown that the electric current density 
carried by an Alfven wing is related to the divergence of the 
electric field through the relation /Yo j ■ V s = S div E where the 
conductance S = (yUoC/i) _1 (l + M 2 A - 2sMa sin#)~'/ 2 , ca being 
the classical Alfven velocity, and Ma = Vq/ca- In this section, 
we derive an analogous relation for the relativistic plasma flow, 
keeping our previous assumptions. From the perfect MHD re- 
lation (fJTJi we calculate div E in the observer's frame Rq'- 



we have 

y 4= ♦(■ 



£ B ■ Vp 

A c 2 

s B ■ v \ Vn s 

7o 



c L c z 



,,2 \ 



^. (48) 
To 



As this vector is constant, its longitudinal and transverse parts 
relatively to the constant vector vo are also two distinct con- 
stant entities. We can recombine the longitudinal and transverse 
components in the following way, 



' sB-v \ s 



,,2 \ 

1 _ J. _ v ° ' Y1 

c 2 c 2 



B 



( Y1) 



+ro 



div E = v • curl B - B • curl v. 
To the first order, 

divE = vo ■ curlBi - Bo ■ curlvi. 

From Eq. d49l , the curl of the velocity is 

s s 
curlvi = — r— curlBj -curl [B (vo • Vj)] 

-LcurlB! - -1-NxBo, 



where 



N = (vo • V)vi + vo x curl vi. 



(56) 
(57) 

(58) 
(59) 

(60) 
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Including Eq. d58b in Eq. 07b . we find 



divE 



vo • 



curlB. 



For further convenience, we note 



U s = v - 



(61) 



(62) 



In Ro, the Alfven wave is stationary, therefore the partial time 
derivatives are null, and Ampere's equation is simply: 



Then, from Eq.doTTi. 



curlB =noj. 



divE =/4)j -Uj 



(63) 



(64) 



Let J s be the projection of the current density along the direc- 
tion of the constant vector U s of equation d49l ). Equation (l64b 
can be written as: 

divE =n J s \U s \. (65) 

Let 6 = n/2 - (Bo, Vo) be the complement of the angle made 
in Ro by the ambient magnetic field and the flow velocity. The 
modulus of Uj (equation d49| >) is 



2s 



u i ~ v o 7 v o B o sin 6» + 



and then: 



J s = divE, 

"/ 2 2 2 \l/2 

1 + C '"V> -iM^sine 



//qCa (l + X 2 - 2sXa sin#) 



1/2 



v *w Vo 2 

i>0 c A 



1 + 



2 2 



2 2 

-^sin, 



(66) 

(67) 
(68) 

(69) 



When ca « c, we are in the conditions studied by Neubauer, 
and we find the same result as in his paper (given in the begin- 
ning of this section). In the case yo >> 1 and ca » c, that is 
relevant for a pulsar's wind, Xa » 1 and 



1 

Hoc' 



(70) 



a conductance which is associated with the impedance of vac- 
uum, equal to H m = hqc — 377 Ohm. The two directions of 
the currents flows correspond to the invariant vectors U + and 
LL.The geometrical configuration of this solution is shown in 
Fig. El 

Now, we ca n adapt quite directly the conclusions of 
Neubaued dl980h to the relativistic inductor. This author con- 
siders a specific model of the wake-aligned currents in the 
Alfven wing, which he assumes to be flowing on the surface of 
an infinite cylinder tangent to the planet's surface, with its axis 
parallel to Alfvenic characteristics. Equation (l67t implies that 
for such currents the divergence of the electric field vanishes 
except on the cylinder's surface. The electric potential can then 
be found by solving Laplace's equation, assuming the electric 



field inside the cylinder to be constant, of intensity This 
simple assumption is motivated by the difficulty to solve for 
the electromagnetic and flow structure in the immediate vicin- 
ity of the solid body. The free parameter £, is the electric field 
along the planet caused by its ionosphere or surface internal 
resistance. At large distances from the wake, the electric field 
converges to the convection field Eo given by the perfect MHD 
relation (fJTJ in the unperturbed wind. Once the electric poten- 
tial is found, the magnetic field and current distribution in the 
wake can be deduced, using in particular equation (|67| >. 

Neubauer gives useful expressions for the total current / 
flowing along an Alfven wing and for the Joule dissipation 
power in the solid body, Ej. Writing Rp for the planet's radius, 
he gets: 

/ = 4(E -E i )RpX A =4l^--E i )R p -Z A , (71) 



Ej = AnR P Ej (Eq - Ej) "La 



(72) 



The Joule dissipation is maximum when = Eq/2. In our es- 
timations, we shall use Neubauer's values for / and Ej. Up to 
unimportant numerical factors, the results (f7Tl -(r72t are simply 
obtained by considering the current / to be driven in a resistive 
load of resistance S^ 1 by a generator of electromotive force 
Uo = 2EqRp applied on two opposite sides of a planet of in- 
ternal resistance %p. Neubauer's parameter E-, is related to Hp 
by: 

El = %p 

E E A ' +K P 



(73) 



It is very difficult to know what the value of Hp really is. 
We therefore regard it, or as unspecified parameters. The 
planet's electrical resistance depends on its constitution, on the 
path of electric currents in it and on the existence or absence 
of some form of ionosphere. The conductivity o f terrestrial sil- 
icate rock is of order cr roci . = 10~ 3 Mho irT 1 (ICookil 1973b . If 
it is assumed that the wake current closes through a layer of 
thickness h at the surface of the planet, the electrical resistance 
Hp of the latter would be of order cr~j Jh, which numerically 
amounts to Hp ~ (/ikm) ■ As soon as h would be larger than 
a few meters, which for this low conductivity is reached in a 
very short time, Hp would be comparable to or less than Hca, 
eventhough the rock intrinsically is a poor conductor. A second 
generation planet, which could be partly metallic, or a partly 
molten body, could have lower resistances. 

6. Discussion and conclusion 

This work shows that a planet orbiting a pulsar develops a sys- 
tem of Alfven wings, caused by its interaction with the sub- 
Alfvenic Poynting-flux-dominated pulsar wind. A system of 
strong electric currents is set. Although this current cannot 
reach the inner pulsar's magnetosphere, it is nevertheless inter- 
esting to compare it to the current at the origin of the pulsar's 
magnetospheric activity, the Goldreich- Julian current Jqj. This 
latter current results from the electromotive field generated by 
the fast rotation of the highly magnetized neutron star and its 
surrounding magnetosphere. For a dipole magnetic field, the 
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Table 1. Input data about the pulsa rs. Rotation period and pulsation, surface magnetic field, mass (in solar masses), radius, and 
Goldreich-Julian current. Source : (ITavlor et al.lEOOOl) catalog, in the SIMBAD database. For PSR 1257+12, the mass and the 
radius are infered from general ideas about pulsars. For PSR 1620-26, the mass has been measured, because the pulsar is in a 
binary system (with a white dwarf companion). 



Name 


P(s) 


a, (s- 1 ) 


B, (Gauss) 


M, (M Q ) 


R, (km) 


Igj (A) 


PSR 1257+12 


0.006 


1010. 


8.8 x 10 8 


1.4 


10. 


4.9 xlO 12 


PSR 1620-26 


0.011 


567 


3. x 10" 


1.35 


10. 


5.4 xlO 12 


PSR 10 ms 


0.010 


628. 


10 8 


1.5 


10. 


2.2 xlO" 


PSR 1 s 


1. 


6.283 


10 12 


1.5 


10. 


2.2 xlO 11 



Table 2. Input data concerning the planets : mass, radius, orbital period, semi-major axis, excentricity. Source : M @ , P or b, a, e 
from "The Extrasolar Planets Encyclopaedia" (http://exoplanet.eu/index.php). The estimate of Rp is based on the hypothesis of 
an average mass density pp = 5000 kg.irT 3 . 



Name 


M P (M e 






Port, (day) 


a(AU) 


e 


PSR 1257+12 a 


0.02 




0.28 


25. 


0.19 





PSR 1257+12 b 


4.3 




1.68 


66. 


0.36 


0.0186 


PSR 1257+12 c 


3.9 




1.62 


98. 


0.46 


0.0252 


PSR 1620-26 a 


794 




9,5 


36367. 


23. 




PSR lsb 10,000km 


3.5 




1.57 


30. 


0.21 


0. 


PSR 10ms b 100 km 


2. x 10" 


-<> 


0.016 


20. 


0.16 


0.3 


PSR 10ms b 1 km 


2. x 10~ 


12 


1.6 x 10- 4 


20. 


0.16 


0.3 


PSR 1 s b 100 km 


2. x 10" 


-6 


0.016 


20. 


0.16 


0.3 


PSR 1 s b 1 km 


2. x 10~ 


12 


1.6 x 10- 4 


20. 


0.16 


0.3 



Table 3. Electric potential drop, total electric current associated to the Alfven wing. Electrical energy Ej max dissipated in the 
Alfven wing. 



Name 


U(V) 


hw (A) 


Ejmax (W) 


PSR 1257+12 a 


1.1 xlO 12 


3.0x10" 


2.5 


xlO 21 


PSR 1257+12 b 


3,5 xlO 12 


9.4 xlO" 


2.5 


xlO 22 


PSR 1257+12 c 


2,6 xlO 12 


7.0 xlO" 


1.4 


xlO 22 


PSR 1620-26 a 


6,0x10" 


1.5 xlO" 


7 


xlO 20 


PSR lsb 10,000km 


3.8 xlO 13 


1.0 10" 


2. 


xlO 24 


PSR 10ms b 100 km 


2.4x10" 


6.X10 6 


1.2 


xlO 16 


PSR 10ms b 1 km 


2.4x1 7 


6.x 10 4 


1.2 


xlO 12 


PSR 1 s b 100 km 


2.4x10" 


6xl0 8 


1.2 


xl0 iu 


PSR 1 s b 1 km 


2.4x10" 


6.x 10 s 


1.2 


xlO 16 



Goldreich-Julian current density is 

J C j = ce £X.B*, 
and the total current is of the order of 



nR PC J C j, 



(74) 



(75) 



where Rpc ~ R?J 2 {£l*l c) 1/2 is the polar cap radius dKirk et al 



2009). Various values of the Goldreich-Julian current Igj are 
given in Table Q] A look at table [3] allows for comparisons 
between Iaw and Iqj. It can be seen that in the case of the 
four known planets, the electric current in the Alfven wings 
is smaller than the Goldreich-Julian current by three orders of 
magnitude. It is also much smaller in the case of small bodies. 
But a planet orbiting a "standard" pulsar with a typical 1 sec- 
ond period and 10 12 G magnetic field would have an Alfven 
wing current I^w of similar amplitude as the Goldreich-Julian 
current Iqj. This is not negligible, when we see that Iqj is the 
basic engine of the pulsar's electrodynamics. 



The practical consequences of such a current are expected 
to be of two kinds : it is shown in a companion paper 



( Mottez & H evvaertsl 1201 lh that it can exert an ortho-radial 
force upon the planet that can, if the magnetic to mechani- 
cal energy coupling is efficient enough, have an incidence on 
the orbit of small circum-pulsar objects such as asteroids or 
comets. The second incidence, more relevant to massive ob- 
jects such as planets, is a possibly associated electromagnetic 
signature, which might be detectable. 

Does the Alfven wing dominates, in terms of energy, the di- 
rect mechanical action exerted by the wind on the companion ? 
To answer this question, we can compare the flux of magnetic 
energy, F m with the flux of mechanical energy, Fk, received the 
companion by direct impingement. The ratio of the Poynting 
flux to the mechanical energy flux is 



Fm 
~F~k 



c (S^M)) 

c (rop c 2 ) ' 



(76) 
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Introducing <x and *P with the Eqs. ( 15161101 ). and \\ )r ~ c, 

Fm _ o"o 
F K To ' 

1 /3 

In the case discussed above where 70 ~ o~ » and ctq >> 1, 



(77) 



Fk 



°? » I- 



(78) 



Therefore, we expect that for a Poynting-flux-dominated wind, 
most of the energy exchange with the pulsar's companion 
comes from the magnetic field. The Eq. (1721 provides a more 
precise insight of what is effectively involved into the wind- 
companion interaction. Considering Ej = Eq and Eq ~ cB^, 



2j}2 



7TR 2 p 



El 2a 



c z B 



Hoc 



Bi 



A'o 



(79) 



This amounts to the totality of the Poynting flux intercepted by 

2/3 

the pulsar's companion. It is therefore larger (by a factor cr Q ' ) 
than the mechanical energy captured by direct impigement. 

This paper rests on the fact that the wind velocity Vo is 
slower than the Alfven wave velocity Va- Otherwise, there 

would be no Alfven wing. In the case of an ideal MHD ra- 

1/3 

dial wind, the Lorentz factor asymptotically approaches cr 
and Eq. (0 shows that the wind remains sub-Alfvenic at any 
distance. If however it is formally considered that the asymp- 
totic Lorentz factor yo<>o scales as cr fl instead, with an expo- 
nent a + 1/3, then if a > 1/3, M Aoo = 1 + cr^ 1 - cr' la ~ 
1 + cr"^ 1 > 1. Therefore, in that case, a transition form a sub- 
Alfvenic to a super-Alfvenic wind occurs at a finite distance. 
Observations of the equatorial sectors of winds driving pul- 
sar wind n ebulae show lower values of the asymptotic mag- 
netization dKennel & Coronitil [l984albl : iGaensler et all 12002) 
and various authors suggest that th e asymptotic value of the 



Lorentz factor is rather y,*, ~ <xo. Arons (2004) argues that 
dissipation must occur in the asymptotic wind zone, in order 
to understand the observed high Lor entz factors and the low 



magnetization. Begelman & Li ( 19941) show that when the flux 



tubes diverge faster than radially, the fast magnetosonic point 
can occur closer to the light cylinder, implying an even closer 
Alfvenic point. In such a circumstance, the existence of Alfven 
wings would depend on the distance from the star to the planet, 
combined with the effect of a non-radial diverging wind flow, 
or dissipation. 

Let us come back to the hypothesis of a planet in a sub- 
Alfvenic wind. The present work provides only orders of mag- 
nitude estimates for the emitted current. The behaviour of the 
Alfven wing at close vicinity of the planet would need a more 
detailed study. We have considered here, as in many other pa- 
pers concerning pulsars, the case of a neutron star magnetic 
field that is aligned with the rotation axis. The study of oblique 
rotators raises more complicated problems. The pulsar wind 
may have different properties, and carry, even in the equatorial 
plane, a non-zero B e component that oscillates at the pulsar ro- 
tation rate Q». An Alfven wavelength (propagating at ~ c) with 
P = 6 ms (case of PSR B 1257+ 12) is of the order of 1800 km, 
which is less than a typical planetary radius. This would put the 
assumption of stationarity into question, though less severely in 



the case of a standard pulsar with P ~ Is, where the wavelength 
(about 300 000 km) would be much larger than the planetary 
radius. 

In spite of the preliminary character of our model, our study 
shows that the consideration of Alfven wings associated to 
planets orbiting pulsars deserve some attention. 

The question of the radio emissions possibly associated to 
the Alfven wings will be addressed in a forthcoming paper. 
Such emissions would provide astronomers with observational 
data relevant to the wind/planet interaction. 
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Fig. 2. The Alfvenic wake of the planet seen from above the 
equatorial plane. 



